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I. INTRODUCTION 

P-fnnctions of a gronp element are the entries of irredncible matrix representations (ir- 
reps) of the element. P-fnnctions of the special unitary group SU(2) are important in 
nuclear, atomic, molecular and optical physics^^-. SU(1,1) is the prototypical non-compact 
semi-simple Lie group, and its R-functions appear in connection with Bogolyubov transfor¬ 
mations, squeezing and parametric downconversion^'^. Methods for construction of intel¬ 
ligent states and the analysis of cylindrical Laguerre-Gauss beams employ R-functions of 
SU(1, P-functions of other Lie groups enable exact solutions to problems in quantum 

optics^ii^. 

One recent application of P-functions of SU(?7,) for arbitrary n is to the BosonSam- 
pling problem, which deals with SU(n) transformations acting on indistinguishable single¬ 
photon pulse inputs^*^. Within the framework of BosonSampling and of multi-photon 
interferometry in general, R-functions provide a deeper understanding of the permutation 
symmetries between the interfering photons. For instance, SU(3) R-functions enable a 
symmetry-based interpretation of the action of a three-channel linear interferometer on 
partially-distinguishable single-photon inputs^*^. Exploiting the permutation symmetries 
present in multi-photon systems reduces the cost of computing interferometer outputs in 
comparison to brute-force techniques^. 

Two existing procedures for computing SU(n) R-function are based on factorization 
and on exponentiation. Both procedures have drawbacks, which we describe as follows. 
Factorization-based methods, which compute SU(n) R-functions in terms of "D-functions of 
subtransformations, are well developed for groups of low rank— . However, generalizing 
these algorithms to higher n requires SU(n — 1) coupling and recoupling coefficients, which 
have limited availability for n > 3, i.e., restricted to certain subgroups of SU(3)^^— . Hence, 
methods for P-functions of higher groups are underdeveloped despite the application of their 
corresponding algebras to diverse problems^I^i^U. 

Another approach for computing SU(?7,) "D-functions involves exponentiating and com¬ 
posing the matrix representations of the algebraiz‘d. This approach has three hurdles. For 
one, this method requires knowledge of all the matrix elements of each generator to be expo¬ 
nentiated. Certain applications require closed-form expressions of "D-functions in terms of 
elements of the fundamental representation; exponentiation-based methods are infeasible for 
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these applications becanse of the difficnlty of exponentiating matrices analytically, especially 
for n > 5. Fnrthermore, if only a limited number of "D-functions are required, exponentiation 
is wasteful because it computes the entire set of "D-functions. 

We overcome the shortcomings of these algorithms by utilizing boson realizations, which 
map the algebra and its carrier space to bosonic operators and spaces respectively. Boson 
realizations arise naturally when considering the groups Sp(2?7,, R), SU(n) and some of their 
subgroups. For instance, SU(1,1), SU(2) and SU(3) boson realizations are used to study 
degeneracies, symmetries and dynamics in quantum systems22,^ii. A wide class of problems 
in theoretical physics rely on boson realizations of the symplectic group^^— . 

Here we aim to devise an algorithm to construct the "D-functions of arbitrary represen¬ 
tations of SU(?7,) for arbitrary n. We approach the problem of limited availability of SU(n) 
D-functions^>^ by presenting (i) a mapping of the weights of an irrep to a graph, (ii) a 
graph-theoretic algorithm to compute boson realizations of the canonical basis states of 
SU(n) for arbitrary n (Algorithm |2] in Subsection IIVBII and (iii) an algorithm that employs 
the constructed boson realizations to compute expressions for "D-functions as polynomials 
in the matrix elements of the dehning representation (Algorithm [3] in Subsection IIV Cj) . 

The rest of the paper is structured as follows. Section [TTl includes dehnitions of the SU(n) 
operators and basis states. In Section IIIIl we dehne SU(n) boson realizations and illustrate 
the calculations of SU(2) "D-functions using SU(2) boson realizations. Section |IV] details our 
algorithms for boson realizations of SU(?7,) basis states and for the SU(?7,) P-functions. We 
discuss potential generalizations of our algorithms in Section |Vl 


II. BACKGROUND: THE SPECIAL UNITARY GROUP AND ITS 
ALGEBRA 

In this section, we recall the relevant properties of special-unitary group SU(n) and its 
algebra su(n). We explain how the 5u{n) D su(n — 1) D ■ • • D su(2) subalgebra chain is 
used to label the basis states of the unitary irreps of SU(n). We present the background for 
n = 2 in Subsection III Al before dealing with SU(n) for arbitrary n in Subsection IIIBl 
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A. SU(2) operators and basis states 


Consider the special nnitary gronp 

SU(2) = {V :V e GL{2, C), CV = 1, det V = 1} 


( 1 ) 


of 2 X 2 special nnitary matrices. Each element of SU(2) can be parametrized by three angles 
= (q!,/ 3,7). The dehning 2x2 representation of an element E(f2) of SU(2) is given by 


e 2 cos I —e sin 


v(f^) = r 1. ^ 

g2d“-7) gi]2 ^ g2^("+7) COS - 


( 2 ) 


2 2 

The Lie algebra corresponding to gronp SU(2) is denoted by su(2) and is spanned by the 
operators Jx, Jy, Jz, which satisfy the angular momentum commutation relations 


[Jxj Jy] '^Jz ) [Jy^ Jz] ^Jx 1 \Jzi Jx\ iJy 
We transform the basis ([3]) of su(2) to the complex combinations 

Jx T '^Jy •} C27 Jx '^Jy 5 J^l 2^/2:, 

which satisfy the commutation relations 


( 3 ) 


( 4 ) 


— 2C'i_2 ; — —2(72,1, [Cl,2, (72,l] — Hi. (5) 

These commutation relations ([5]) facilitate the construction of a (2J + l)-dimensional irrep 
with carrier space spanned by basis states {\J,M) : —J < M < J}^. The integer 2M is 
the weight of the eigenstate \J,M) for 


Hi \J,M) = 2M|J,M). 


( 6 ) 


The operators (7i,2 and (72,i act on eigenstates of Hi by raising or lowering the weight 2M 
of the states 


(7i, 2|J,M) = v/j(jTiy^WWXMTl) |J,M + 1), (7) 

(72,1 1^, M) = ^/J{J + 1) - M(M - 1) I J, M - 1), (8) 

where 2J is the highest eigenvalue of TTi. 

Each basis state of a hnite-dimensional irrep of SU(2) is labelled by integral weight 2M G 
{—2J, —2J + 2,..., 2J — 2, 2J}. The unique basis state | J, J) is called the highest-weight 
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11 _i\ 11 i\ 

I 2’ 2/ I 2’ 2/ 

•- • 

Hi |2,-2) |2,-1> |2,0> |2,1> |2,2) 

^2,1 i -•-> Cl^2 •-•-•-•- • 

(a) (b) 

FIG. 1. (a) Generators of the su(2) Algebra. The action of the raising and lowering operators 
Cl, 2 ; (^ 2,1 on the basis states is represented by the directed lines. The basis states are invariant 
under the action of the Cartan operator Hi, which is represented by the dot at the centre, (b) 
SU(2) irreps labelled by highest weights 2M = 1 and 2M = 4 respectively. The dots represent the 
basis states while the lines connecting the dots represent the transformation from one basis state 
to another by the action of the su(2) raising and lowering operators. The red dot represents the 
hws, which is annihilated by the action of the raising operator Ci^ 2 - 

state (hws) and is annihilated by the action of the raising operator C'i_ 2 . The representation 
is labelled by the largest eigenvalue 2J of Hi. Basis states of an SU(2) irrep are visualized 
as collections of points on a line with the location of each point related to the weight of 
the state. Figure [U gives a geometrical representation of the action of su(2) operators and 
illustrative examples of SU(2) irreps. 

B. Basis states and P-functions of SU(n) for arbitrary n 

Next we consider the case of arbitrary n. The unitary group U(n) is the Lie group of 
n X n unitary matrices 


U(n) := {I/: 1/e GL(n,C),I/V = 1}. (9) 

The corresponding Lie algebra is denoted by u(n). The complex extension of u(n) is spanned 
by operators {Cij: i,j G 1, 2,... n} satisfying the canonical commutation relations 

[Cjj, Ck,i] = Sj^kCi^i — di^iCkj. (10) 

The group SU(n) is the subgroup of those U(n) transformations that satisfy the additional 
property det V = 1; i.e., 


SU(n) := {V: V e U(n),detl/ = 1}. 


( 11 ) 
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The U(n) "D-fnnctions differ from the SU(n) "D-fnnctions by at most a phase, and we con¬ 
centrate here on the SU(n) case. 

The operator N = Ci^i + 6 * 2,2 + • • • + is in the centre^ of u(n). The Lie algebra su(n) 
is obtained from u(n) by eliminating the operator N . The n — 1 operators 

H, = a,, - Q+ 1,,+1 V* e { 1 , 2 ,..., n - 1 } ( 12 ) 

commnte with each other and span the Cartan snbalgebra of su(n). Hence, we have the 
following dehnition of the su(n) algebra. 

Definition 1 (su(n) algebra^). The algebra 5u{n) is the span of the operators {Cij\ i,j G 
{1, 2,..., n}, i ^ j} and {Hi : Hi = Ci^i — i G {1, 2,..., n — 1}} where the operators 

{Cij} obey the commutation relations 

[Ci,j, Ckf = dj^kCi^i — di^iCkj- (13) 


The linearly independent (LI) su(n) basis states span the carrier space of su(n) representa¬ 
tions. Each basis state is associated with a weight, which is the set of integral eigenvalnes 
of the Cartan operators. 


Definition 2 (Weight of su(n) basis states^). The weight of a basis state is the set A = 
(Ai, A 2 ,..., A„_i) of n — 1 integral eigenvalues of the Cartan operators [Hi, H 2 ,..., 
su(n) basis states have well defined weights. 


Of the — 1 elements, n — 1 Cartan operators generate the maximal Abelian snbalgebra 
of su(? 7 ,). The remaining operators satisfy the commntation relation 




{ Pi,jkCj^kt < /c. 


(14) 


for Cartan operators Hi of Dehnition [T] and for positive integral roots (dijk- The opera¬ 
tors {Cj^k'- j < k} dehne a set of raising operators. The remaining off-diagonal operators 
{Cj^k '■ j > k} are the su(n) lowering operators. Each irrep contains a nniqne state that has 
nonnegative integral weights K = (ki, ..., Kn-i) and is annihilated by all raising operators. 
This state is the hws of the irrep. 
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Definition 3 (Highest-weight state). The hws of an SU{n) irrep is the unique state that is 
annihilated according to 


Ci.i I'lpL,} =0 Vi < j. !,J e {1.2, ....n} 


(15) 


by the action of all the raising operators. 

The weight of the hws also labels the irrep; i.e., two irreps with the same highest weight 
are eqnivalent and two eqnivalent representations have the same highest weight. Hence, we 
label an irrep by A' = (ki, K 2 , , Kn-i) if the hws of the irrep has weight A = K. 

Whereas in SU(2) the weight 2M and the representation label J are enough to uniquely 
identify a state in the representation, this is not so for SU(n) representations. In general, 
more than one SU(?7,) basis state of an irrep could share the same weight. For example, 
certain states of the K = (2,2) irrep of SU(3) irrep have the same weight (Fig. [2]). The 
number of basis states that share the same SU(n) weight A = (Ai, A 2 ,..., A„_i) is the 
multiplicity M(A) of the weight^. Hence, uniquely labelling the SU(n) basis states requires 
a scheme to lift the possible degeneracy of weights. 

One approach to labelling the SU(n) basis states involves specifying the transformation 
properties under the action of the subalgebras of su(n). We restrict our attention to the 
canonical subalgebra chain 


5Ui,2,...,n(n) D SUi,2,...,n-l(n - 1) D • ' ' D SUi,2(2), 


(16) 


where 5Ui^2,...,m(^) is the subalgebra generated by the operators {Qj: i, j G {1,2,..., m} , i ^ 
j} and {Hk \ k G {1,2,... ,m — 1}}. Details about the choice of subalgebra chain are pre¬ 
sented in[^ Henceforth, we drop the subscript and denote sui, 2 ,...,m("i) by su(m). 

The canonical basis comprises the eigenstates of the 5u{m) generators for all m < n 
according to the following dehnition. 

Definition 4. (Canonical basis states) The canonical basis states of SU{n) irrep are 


those states 



( 17 ) 


that have well defined values of 


1. irrep labels for su{m) algebras for all {m : 2 < m < n} and 
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(^ 2,3 C'l.S 



FIG. 2. (a) Generators of the su(3) algebra. The action of the raising operators {C'i^2j C'l.s, <^2,3} 
and lowering operators {C'2,1, 6*2,3} on the canonical basis states and their linear combinations 

is represented by the directed lines, (b) The SU(3) irrep labelled by its highest weight (ki,K 2) = 
(2,2). The dots and circles represent the canonical basis states. The dimension of the space of 
states at a given vertex is the sum of the number of dots and the number of circles at the vertex, 
for instance weights associated with dimension two are represented by one dot and one circle. The 
lines connecting the dots represent the transformation from states of one weight to those of another 
by the action of SU(3) raising and lowering operators. The red dot represents the highest weight 
of the irrep. A unique hws occupying this weight is annihilated by the action of each of the raising 
operator. 


2. su(m) weights i.e., eigenvalues of the Cartan operators of su{m) algebras for all 
{m : 2 < m < n}. 

Consider the example of the (ki, ^ 2) = (1,1) irrep of SU(3). There are two basis states with 
the weight (Ai, A2) = (0, 0). We can identify these two states by specifying 

1. the su(3) irrep label = (ki, K2) = (1,1) and the su(2) irrep label = (ki) = (0) 
or = (ki) = (1 ). 

2. the su(3) weights A*^^) = (Ai, A2) = (0,0) and su(2) weight A*^^^ = (Ai) = (0). 
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The connection between onr labelling of canonical basis states of Dehnition 0] and the 
Gelfand-Tsetlin patterns^ is detailed in[Bl The canonical basis state for 

which for all m G {2,..., n} is the highest weight of the irrep 

The relative phases between the canonical basis states are hxed by comparing with the 
phase of the hws^. Matrix elements of the simple raising operators i G {1,..., n — 1} 

are set as positive^. Thns, we impose the following additional constraint on the canonical 
basis states 


/./, ^1.2P2,3.. P: 

y V^hws 2 02,3 


'n — l,n 

n—l,n 


I K(^\...,K(^\kW\ ^ PI 

V’aM A(2) / — 


(18) 


for all canonical basis states, for positive integers pe^e+i- 

P-fnnctions are the matrix elements of SU(n) irreps. The rows and columns of SU(?7,) 
matrix representations are labelled by SU(n) basis states. The expression for SU(n) "D- 
functions generalize those of the SU(2) D-functions (l25j) with M, M' replaced by suitable 
labels for weights and J replaced by suitable subalgebra labels. 


Definition 5 ("D-functions). V-functions of an SU{n) transformation 1/(12) are 

^A(") ,...,A(3), A(2) ;A'(") ,...,A'(3), A'(2) ■“ \^A(") ,...,A(®/A(2) 

where 12 = {ci;i,ci;2 ,... ,a;„2_i} is the set of — 1 independent angles that parameterize an 
SU{n) transformatior^. 


”^A'(") ,...,A'(3), A'(2) 


(19) 


Note that SU(n) 12-functions (IT^ are non-zero only if the left and the right states belong 
to the same SU(n) irrep, i.e.. 


J^{n) ^ 


.t-,A'("),...,A'( 3) pp(2).jp/(n)_ pp/(3) pp/(2) _ 

^A(") ,...,A(3), A(2) ;A'(") ,...,A'(3), A'(2) ^ 


( 20 ) 


I2-functions of an irrep K refer to those "D-functions for which A'G) = 

We approach the task of constructing SU(n) "D-functions by using boson realizations of 
SU(?7,) states. In the next section, we dehne boson realizations and illustrate the construction 
of SU(2) 12-functions using SU(2) boson realizations. 


III. BACKGROUND: BOSON REALIZATIONS OF SU(n) 

In this section, we describe boson realizations, which map su(n) operators and carrier- 
space states to operators and states of a system of n — 1 species of bosons on n sites 
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respectively. We first present the mapping for n = 2 and illnstrate SU(2) D-fnnctions 
calculation using the SU(2) boson realization in Subsection IIIIAI Boson realizations of SU(n) 
for arbitrary n are dehned is Subsection IIIIB1 


A. SU(2) boson realizations 


The commutation relations ([5]) of {Ci 2, C2,i, hfi} are reproduced by number-preserving 
bilinear products of creation and annihilation operators that act on a two-site bosonic system. 
Specihcally, the su(2) operators have the boson realization 

C*i,2 Ci^2 ■= ci\(i2 , C*2,i ^ C2,i := a^tti, Hi 1 -^ hi ■= a\ai — 0 . 2 ^ 2 ; (21) 


where the bosonic creation and annihilation operators obey the commutation relations 


ai,a] 




ala] 


= 0 . 


( 22 ) 


Here and henceforth, we use lower-case symbols for boson realizations of the respective 
upper-case symbols. Explicitly, 


[hi, 01^2 ] — 2Ci^ 2 [hi, C2,l] — —2C2,1 [Ci, 2; C2,l] — hi. 


(23) 


The operators {ci,2, C2,i, hi} also span the complex extension of the su(2) Lie algebra. 

Boson realizations map the states in the carrier space of SU(2) to the states of a two-site 
bosonic system. Specihcally, each basis state of the (2J -|- l)-dimensional SU(2) irrep maps 


\J,M) ^ 


(a 


Uj+MfJ\J-M 


1 ) 


do 


yJ{J + M)\{J - M)\ 


| 0 ) 


(24) 


to the state of a two-site system with J + M and J — M bosons in the two sites respectively. 

The (2 J -|- l)-dimensional irreps of SU(2) map to number-preserving transformations on a 
two-site system of 2J bosons in the basis of Eq. ([24)). The elements of these (2J-I-1) x (2J-I-1) 
matrices are the SU(2) "D-functions 


■■= {J,M'\V(Sl)\J,M) 


(25) 


for irrep J and row and column indices M', M. The expression for "D-functions fl2^ of SU(2) 
element V (fl) can be calculated by noting that the creation operators transform under the 
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action of V of Eq. (j2]) according to 


a\ —y Viiaj + Vi2Qj\, 

a\ V2iCi| + V22Ci\, (26) 


where V is the 2x2 fnndamental representation of The state \J,M) fl2^ thus 

transforms to 


\J,M) 


N J+M / \ J-M 

VllCLl + Vi2n2 ) ( T ^22^2 


| 0 ) 


(27) 


^{J + M)\{J - M)\ 

as the vacuum state |0) is invariant under the action V. Using Eqs. fj2Tll and fl271) . we obtain 


D 


j 

M'M 


(fi) = 0 


a 


J+M'J-M' 


1 


dr) 


xJ+Af/ .\J-M 

Uiiaj + Vi2a2 J + U22a2 


^(J + M')!(J - M')\^J{J + M)!(J - M)\ 


(28) 


which can be evaluated using the commutation relations of the creation and annihilation 
operators ([22]).— 

In this paper, our objective is to generalize Eqs. (12T]) and fl24l) systematically from n = 
2 to arbitrary n. In the next subsection, we dehne boson realizations of operators and 
carrier-space states of su(?7,). Furthermore, we construct the boson realization for the hws 
of arbitrary SU(n) irreps. 


B. SU(n) boson realizations for arbitrary n 


SU(n) boson realizations map SU(n) states and su(n) operators to 

states and operators of a system of bosons on n sites. Bosons are labelled based on the site 
i G {1,2,... ,n} at which they are situated and by an internal degree of freedom, which 
is denoted by an additional subscript on the bosonic operators. The bosonic creation and 
annihilation operators on this system are 


|a|j: i G {1,2,... ,77,}, j G {1,2,... ,n - 1}| Creation (29) 

: /c G {1, 2,..., n}, / G {1, 2,..., n — 1}| Annihilation, (30) 


where the hrst label in the subscript is the usual index of the site occupied by the boson. 
The second index refers to the internal degrees of freedom of the boson. Each boson can 
have at most n — 1 possible internal states to ensure that basis states can be constructed for 
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arbitrary irreps. In photonic experiments, this internal degree of freedom conld correspond to 
the polarization, freqnency, orbital angular momementum or the time of arrival of photons. 

The su(n) operators are mapped to number-preserving bilinear products of boson creation 
and annihilation operators. Specihcally, raising and lowering operators Cij of su(n) map to 
bosonic operators according to 

n—1 

Cij Cij ■= ^ ^ ^\,k^j,k- (31) 

k=l 

Operators {cij} make bosons hop from site j to site i. The operators hi are the image of 
the Cartan operators Hi: 

Hi I ^ hi — Oj^di (32) 

Operators {hi} count the difference in the total number of bosons at two sites and commute 
among themselves. As usual, we used the upper-case symbols to denote the su(n) elements 
and the corresponding lower-case symbols for the respective boson operators. 

The boson realizations of the basis states of SU(n) are obtained by the action of poly¬ 
nomials in creation operators {ajj : i G {1, 2,..., n},j G {1, 2,..., n — 1}} on the n-site 
vacuum state |0). Each term in the polynomial is a product of 

Nk = Ki 2K2 • -f- (n — l)Kn—l (33) 


boson creation operators for basis states in irreps K = {ki, K 2 , ■ ■ ■, Hn-i)- Therefore, an 
SU(n) basis state is specihed by the coefficient of a polynomial consisting of terms that are 
products of Nk creation operators. 

The hws of a given SU(n) irrep can be explicitly constructed in the boson realization (as 
polynomials in creation and annihilation operators) according to the following lemma. 


Lemma 6 (Boson realization of hwa^*^). The bosonic state 


( .t 


\i^Ls) = det 


1,1 


\ 


Hn — 1 


■■■det ( “j’M det(^a^) |0) (34) 


1,1 • • • ®n—l,n—1 J 

is a hws for a given SU{n) irrep K = (ki, K 2 , ■ ■ ■, n,n-i)- 
One can verify that the state |V’^g) fl3T)) is annihilated 

Cj,k =0 Vj < A; 


®2,1 ®2,2 


(36) 
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by the action of any of the raising operators. 

Thns, the hws of any irrep can be constrncted analytically nsing Lemma [6l In the 
following section, we provide an algorithm to constrnct each of the basis states of arbitrary 
SU(?7,) irreps. Fnrthermore, we present an algorithm to compnte expressions for SU(n) P- 
fnnctions in terms of the entries of the fnndamental representation. 

IV. RESULTS: ALGORITHMS FOR BOSON REALIZATIONS OF SU(n) 
STATES AND FOR SU(n) P-FUNCTIONS 

In this section, we present three algorithms^. Algorithm [1] (boson-set algorithm) con- 
strncts basis sets for each weight of a given su(n) irrep. Algorithmic] (canonical-basis-states 
algorithm) employs the boson-set algorithm to compnte expressions for the canonical basis 
states of a given SU(n) irrep K. The states thus constructed are used by Algorithm [3] to 
calculate the "D-functions of a given SU(n) transformation. 

Algorithms [T] and 2 rely on mapping the SU(n) irrep to a graph and systematically 
traversing the graph to obtain basis states. The vertices of the irrep graph are identihed 
with the weights of the given irrep of SU(n) and the edges with the action of the elements 
of the Lie algebra su(n) on the states. Specihcally, the irrep graph G = {y,£) of an SU(n) 
irrep is dehned as follows. 

Definition 7 (Irrep graph). The bijection 

V: {Ai, A 2 ,..., Arf} V (36) 

maps the set {Ai, A 2 ,..., A^} of the d weights in the given irrep to the vertices 

V = {v{Ai),v{A2),...,v{Ad)} (37) 

of its irrep graph. Vertices v(Ak) andv{Ai) are connected by an edge Cj = (n(Afc), ^(A^)) G £ 
iff 3 Cij,Ak,Ai such that 

= \'fA,) , (38) 

where It/'Afe) and |V’a^) are SU(n) states that have weights Ak and A^ respectively. In general, 
states iV'Afc) and It/^A^) are linear combinations of canonical basis states. Edges £ together 
with the vertices V define the irrep graph G = {y,£). 
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FIG. 3. The first step of the basis-set computation algorithm (Algorithm [T|) illustrated for the 
(2,2) irrep of SU(3), where the dimension of the space of states at a given vertex is the sum of 
the number of dots and the number of circles at the vertex. The algorithm constructs the hws 
(occupying the red vertex) using Lemma [6l The lowering operators can transform states at one 
vertex to states at another vertex along different paths connecting the starting and the target 
vertex, for instance the two paths coloured green and blue. Lowering along the different paths 
to reach a target vertex will generate the same number of LI as the weight multiplicity. In our 
illustration, we obtain a basis set that contains two independent states at the target vertex. The 
algorithm traverses the irrep graph systematically until all basis sets are calculated. 


More than one basis state can have the same weight. The nnmber of basis states sharing 
a weight Aj is defined as the mnltiplicity M(Aj) of the weight. In other words, each vertex 
n(Aj) is identihed with an M(Aj)-dimensional space spanned by those canonical basis states 
that have weight Aj. The vertex space and vertex basis sets are defined as follows. 


Definition 8 (Vertex spaces). The vertex space of v{Ai) is the span 


^(Aj) = span (pj.) , , 



(39) 


of the canonical basis states (Definition^ that have the weight Ki. 



spans the vertex space \l/(Aj) of n(Aj). 
follows. 


of canonical basis states is not the only set that 
In general, basis sets of T(Aj) can be dehned as 
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Definition 9 (Vertex basis sets). The set 






dMiAi)) 


(40) 


is called the basis set of a vertex v{Ai) if it spans the vertex space tt'(Ai) ofv{Ai), i.e., 






= ^(A. 


(41) 


The states 


*a’/. 




0^(A“)) ) [> are linear combinations of the canonical basis 


states 


, 4’Ai'^ ) • • •) Algorithm [T] CO mpntes basis sets of the spaces 4/(Ai 


for each of the d weights Aj that occurs in a given irrep. 


A. Basis-set algorithm (Algorithm [T]) 


The basis-set algorithm, which hnds the basis sets for a given SU(m) irrep, is the key 
subroutine of our canonical-basis-state algorithm. Algorithm [1] requires inputs 
m, where IV'ws) ^ irrep K of su(m) algebra. The state IV'w) ^ bosonic state, 

which is expressed as a summation over products of 115^ creation operators {a|^: i G 
{1, 2,..., m},j G {l,2,...,m — 1}}. This summation acts on the m-site bosonic vacuum 
state to give an A^-boson state. The algorithm returns multiple sets 







IV')} 


,(«(».)))} 


5 







,,«(A.„)} 


(42) 


of su(m) states, with each set spanning the space \I/(Aj) (15^ at a different vertex n(Aj) 
in the SU(m) irrep K. The states in the output basis sets are represented as polynomials 
in lowering operators acting on the hws, or equivalently as polynomials in creation and 
annihilation operators acting on the n-site vacuum state. Figure [3] is an illustrative example 
of the algorithm. 

A modihed breadth-hrst search (BFS) graph algorithm^^— is used to traverse the irrep 
graph for states. As in usual BFS, we maintain a queued, called currentQueue, of the states 
that have been constructed but whose neighbourhood is yet to be explored. The algorithm 
starts with the given hws in currentQueue and iteratively dequeues a state from the front 
of the queue. States neighbouring the dequeued state are obtained by enacting one-by-one 
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Algorithm 1 Basis-Set Algorithm 

Input: 

• hws IV'ws) ^ Degree Nk (l33l) polynomial in bosonic creation operators. 

• m € Z+ t> is a hws of SU(m) irrep K. 

Output: 

• {Ai, A 2 ,... ,Ad: Ai e (Z+ U 0)"*“^} > List of weights in the irrep graph of K. 

• d, Basis sets (H2|) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 


{|< 
(lAAlA’) 


a; 

‘S' 


Am(a,)) 

^Ai 


}’{ks 




i(M(*2))\ 




/ 




‘S 


‘S) 






}■ 


hws) 


procedure BASisSET(m, 

Initialize empty statesList, empty weightList and currentStateQueue •(— |'k’] 
while cnrrentStateQuene is not empty do 

cnrrentState ^ DEQUEUE(cnrrentStateQuene) 
for CurrentOperator € set of su(m) lowering operations do 
newState ■(— CuRRENTOPERATOR(cnrrentState) 
if newState 7^ 0 then 

if weight of cnrrentState is already in stateList then 

if cnrrentState is LI of stateList states with same weight then 
independentState ^ NORMALiZE(newState) 

Enqueue independentState in currentStateQueue 
Add independentState to stateList 
Add weight of independentState to weightList 
end if > Else, do nothing, 

else 

Enqneue newState in cnrrentStateQuene 
Add {weight(newState),newState} to stateList 

end if 
end if 
end for 
end while 
Retnrn stateList 
end procedure 


each of the lowering operators of the algebra. The newly found states are enqueued into the 
rear of currentQueue, and the current state and its weight are stored. 

We modify BPS to handle vertices with weight multiplicity greater than unity as follows. 
While traversing the irrep graph, the algorithm directly enqueues the hrst state that is found 
at each vertex. When the same vertex is explored along a different edge, i.e., by enacting 
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different lowering operators, a different state is fonnd in general. If the newly constrncted 
state is LI of the states already constrncted at the vertex, then the new state is enqneued 
into currentQneue. 

The algorithm trnncates when a state in currentQneue is annihilated by all of the lowering 
operators and there is no other state in the queue. This hnal state must exist because the 
number of LI states in a given SU(n) irrep is hnite according to the following standard result 
in representation theory. 

Lemma 10 (Dimension of an SU(n) irrep2^). The dimension of the carrier space of an 
SU{n) irrep K is 

Ak ^=M(Ai) + M(A2) + ■ • ■ + M(Arf) 

= (1 + K,) (1 + K.) . . . (1 + K„_,) (l + (l + 

"A 2 )[ 3 )[ 3 ) 

. l^n-3 + ^n-2 + ^n-1 \ A + ^2 + ' ' ‘ + ^n-1 

3 / V n-1 

Now we prove that the basis-set algorithm terminates. The proof relies on the fact 
that the carrier space of SU(m) irrep is hnite-dimensional 1 Lemma [T0|) . The algorithm’s 
computational cost is quantihed by the number of times the lowering operators are applied 
on the hws or on states reached by lowering from the hws. We show that the computational 
cost of Algorithm [T] is linear in the dimension Ak of the irrep whose hws is given as input 
and polynomial in n. 

Theorem 11 (Algor ithm[T] terminates). Suppose Algorithm{l\ receives as input an hws 
of an SU{m) irrep K. Then the algorithm terminates after no more than AKm{m — l)/2 
applications of lowering operators. 

Proof. The proof is in two parts. Firstly, the number of states that enters currentStateQueue 
is bounded above by the dimension Ak (|^ of the irrep space. Secondly, as each state that 
enters currentStateQueue is acted upon by no more than n{n — l)/2 lowering operators, the 
number of lowering operations performed is less than or equal to Axn(n — l)/2. 

We show that the number of states that enter currentStateQueue is no more than Ak as 
follows. As each currentState is a linear combination of states obtained by acting lowering 
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operators (Line [6]) on the given hws, each state that enters cnrrentStateQnene is in the irrep 
labelled by the hws. Moreover, each state entering the qneue is tested for linear independence 
(Line [9]) with respect to the states already obtained. Any state that is not LI is discarded. 
Therefore, each enqnened state (Algorithm HI Line [T6|l is in the correct irrep K and is LI of 
each other enqnened state. Thus, the number of states that ever enter currentStateQueue 
is no more than the number of LI states in irrep K. 

In each iteration of the algorithm, we act all the lowering operators on the states in cur¬ 
rentStateQueue. The number of lowering operations is thus bounded above by the product 
AKn{n — l)/2 of the number of states that enter currentStateQueue and of the number of 
lowering operators in the su(n) algebra. The algorithm thus terminates after no more than 
AKn{n — l)/2 applications of lowering operators. ■ 

We now prove that the algorithm returns the correct output on termination. The proof 
requires the following lemma stating that each canonical basis state can be obtained by 
enacting only with the lowering operators on the hws. 

Lemma 12 (Every basis-state can be reached by lowering from the hws^). No canonical 
basis-state of a given SU{n) irrep K is LI of the states obtained by lowering from the hws 
by the action 

^ikJk ■ ■ ■ l^/^hws) ii < i < k (44) 

of k < number of 5u{n) lowering operators on the hws of the irrep. 

Lemma [I2] implies that each basis state can be constructed by linearly combining states 
obtained on lowering from the hws. Algorithm [T] leverages from the construction of Eq. fl44|) 
and from testing linear independence to construct the basis sets. 

The correctness of the basis-set algorithm is proved as follows. We show that each state 
obtained by enacting any number of lowering operators on the hws is LD on the states 
returned by the algorithm. Each canonical basis state is LD on the states obtained by 
lowering from the hws in turn, so each canonical basis state is LD on the algorithm output. 
The algorithm only constructs states in the correct irrep so Algorithm [T] returns a complete 
basis set at each weight of the irrep on truncation. 
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Theorem 13 (Algorithm [T] is correct). The sets 



of states returned by Algorithm\l\ span the respective vertex spaces \t^(Aj) / T^) at each ver¬ 
tex Aj of the given irrep K. 

Proof. We hrst prove by indnction that each state in the form of Eq. 0441) is LD on states 
in the algorithm ontpnt. Onr indnction hypothesis is that each state 

^ikdk ' ' ' IV^hws) ) (45) 

which is obtained by acting i lowering operators on the hws, is LD on the states retnrned by 
the algorithm G The proof of the hypothesis follows from mathematical indnction 
over i. 

The indnction hypothesis is trne for base case k = 1. In the hrst iteration, the algorithm 
enacts all the lowering operators on the hws (Algorithm [1] Line |6]) and saves each of the 
obtained states. No /c = 1 state fH5l) is omitted because the vertices neighbouring the hws 
vertex are all being explored for the hrst time. Hence, all the states that can be reached by 
lowering once from the hws are added to currentStateQueue and, eventually, to stateList. 

Assume that the induction hypothesis holds for k = i, i.e., each k = £ state is LD on 
the states in stateList. We prove that the hypothesis holds ioi k = £ + Ihj contradiction. 
Suppose there exists a state that can be reached by enacting £ + 1 lowering operators on the 
hws but is LI of stateList. Let l"^) = • • •CiajaDiji It/^hws) be such a state. 

Consider now the state \ip) = • • •Cj2j2Di,ji iV'hws) obtained by enacting one less low¬ 
ering operation from the hws; i.e., I"^) = Iv^)- We have assumed that the induction 

hypothesis holds for k = £. Therefore, \ip) is LD on the states constructed by a algorithm. 
In other words, 

= ( 46 ) 

is LD on the stateList elements {|0j) : j G {1, 2,..., J}} for complex numbers Oj. 

The algorithm enacts the lowering operator on each \(j)j) and the resulting states 

are either stored in stateList or are LD on elements in stateList. Therefore, the elements of 
the set \(j)j) : j G {1,2,..., J}} are LD on the elements of stateList. Hence, the 
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element \^p) is also LD on the elements of stateList. This dependence contradicts 

the snpposition that IV') = \^) is LI of stateList, thereby proving the induction 

hypothesis for k = ^ + 1. 

The induction hypothesis is true for £ = 1 and is shown to hold for A; = £ + 1 if it holds 
for k = £. Thus, our induction hypothesis is true for all i G Z’*'. Every state obtained of 
irrep K obtained by lowering from the hws is linearly dependent (LD) on the basis sets that 
are returned by the algorithm. 

We know from Lemma [l2] that each canonical basis state is LD on the states obtained 
by lowering. Hence, each canonical basis state is LD on the states obtained at the output 
of the algorithm. Therefore, the state returned by the algorithm span the space of irrep K 
states, and the output basis sets span the set of all basis states of the given irrep K. ■ 

We have proved that Algorithm [T] terminates and that it returns the correct basis sets 
on termination. Now we present our algorithm for the construction of the canonical basis 
states. Furthermore, we prove the correctness and termination of the canonical-basis-states 
algorithm. 


B. Canonical-basis-states algorithm (Algorithm [2]) 

The algorithm for constructing the canonical basis-states of SU(n) requires inputs n G Z+ 
and the irrep label K. The algorithm returns expressions for all the canonical basis states 
in the given irrep. Figure 0] illustrates SU(3) basis-state construction using our algorithm. 
Algorithm [2] details the step-by-step construction of the canonical basis states. 

The canonical-basis-states algorithm proceeds by partitioning su(n) basis sets into su(m) 
basis sets for progressively smaller m over n — 1 stages. In the hrst stage, the algorithm 
employs Lemma [6] to construct the hws of the given irrep K (Algorithm [2l Line [2]). Algo¬ 
rithm [T] is then used to construct the basis sets of the SU(n) irrep of the constructed hws 
(Line [3]). 

By the (n — m)-th stage, the algorithm has partitioned the entire su(n) space into basis 
states of the SU(m -I- 1) irreps. In this stage, each of the su(m -|- 1) basis sets is partitioned 
into su(m) basis sets by using 5u(m) operators. The algorithm searches each SU(m -|- 1) 
irrep graph for the vertex that has the highest multiplicity. 
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FIG. 4. Diagrammatic representation of the main algorithm for n = 3. The dots and circles 
represent the canonical basis states. The dimension of the space of states at a given vertex is the 
sum of the number of dots and the number of circles at the vertex, for instance weights associated 
with dimension two are represented by one dot and one circle. The lines connecting the dots 
represent the transformation from states of one weight to those of another by the action of SU(3) 
raising and lowering operators. We use Algorithm [1] to construct basis sets for each vertex in the 
SU(n) irrep graph. Once the basis sets for the SU(n) irreps are computed, the algorithm enacts the 
su(n — 1) raising operators on the (n — l)-dimensional sub-irreps to find the su(n — 1) hws. Then 
the algorithm starts with the 5u(n — 1) hws and employs the basis-set construction (Algorithm 1) to 
find all the states in the su(n — 1) irrep labelled by the hws. The states thus obtained are subtracted 
from the set of su(n) states. A new state is chosen from the weight of highest multiplicity and the 
process repeated until all the su(n — 1) irreps are found. 

An arbitrary linear combination of the basis states at this vertex is chosen. The algorithm 
then enacts all the raising operators in the su{m) subalgebra on this linear combination until 
the action of each of the raising operators annihilates the state. The state thus obtained 
is the hws of an SU(m) irrep, whose label can be calculated by enacting the Cartan 
operators on the state. 

Next the algorithm performs the basis-set construction algorithm on the su(m) hws em¬ 
ploying only the su(m) lowering operators. This procedure gives us sets of basis states that 
belong to the SU(m) irrep The irrep basis sets are stored and are then sub¬ 

tracted from the SU(m -|- 1) states. The algorithm iteratively (i) starts from the highest 
multiplicity vertex of SU(m-l-l) irrep graphs, (ii) constructs a hws by raising, (iii) stores the 
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Algorithm 2 Canonical-basis-states algorithm 


Input: 


n 


G Z+ 


K = , Hn-i) G u {0}) 


> Algorithm constructs basis sets of su(n) algebra. 

i> Label of SU(n) irrep. 


71—1 


Output: 


{ ( ’ '^(3^'^( 2 / ) ’ • • • 7 ..., I t> List of all canonical basis states 

and weight labels in the irrep = K. 


9: 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 

22 

23 

24 

25 


procedure CANONiCALBASisSTATES(n, K) 

Initialize empty basisStatesList, hws •<— 

SUmStates, SUnStates •(— BASisSET(n,hws) 
while SUnStates is not empty do 
for m G {n, n — 1,..., 2} do 

Amax su(m) weight with highest number of states in SUmStates. 

V’r^x^ arbitrary superposition of states at Amax in SUmStates. 

Apply su(m — 1) raising operators on V’r^x^i reach su(m — 1) hws 'i/’hws 
i)^Weight(|4™-'))). 

SUmStates ^ BASisSET^m — 1, V’hws"^^))- 

if m = 2 then 

for All states \'ip) in SUmStates do 

{Ai”i, ..., A(^)} ^ Weights(IV')) 0 su(m) weights Vm < n . 

Concatenate (lip) ] K^'^\ ..., A®; Ai""),..., A^^i) basisStatesList 
Subtract SUmStates from SUnStates. 

end for 

end if i> Else, do nothing. 

end for 
end while 

for All states in statelist do 


Act {Cl, 2 , C2,3, ■ ■ ■, Cn-I,n} on until hws 




for 


I^W) 

end for 

Return basisStatesList 

end procedure 


Ip. 


(h \_ 

hws /~ 


''/'hws) reached. 


kws>- 


basis sets of SU(m) irreps corresponding to this hws and (iv) subtracts them from SU(m-l-l) 
states until all the states in the su{m + 1) are partitioned. 

At the end of n — 1 stages, we have a list of basis sets of the SU(n — 1) irreps. We 
iteratively perform the process of hnding basis sets for smaller subgroups until we reach 
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SU(2) basis sets, which are known to have nnit mnltiplicity. Hence, the algorithm retnrns 
the basis states that are eigenvectors of the Cartan operators of all SU(m) : m < n groups. 

The relative phases between the basis states are hxed by imposing Eq. flTSll . Each of the 
constructed basis states is acted upon by the simple raising operators {Ci, 2 , 1 ^ 2 , 3 , • • •, C'n-i,n} 
until the hws is reached. The phase of this hws obtained by raising is required to be the 
same for all basis states. Our algorithm multiplies each of the basis states by a phase factor 
(Line [22]) to impose the phase convention Eq. flT8|) and returns the set of canonical basis 
states. 

Now we prove that the canonical basis states algorithm terminates. The proof of termi¬ 
nation uses the facts that the number of basis states is equal to the dimension Ak of the 
irrep and that each basis state is added to currentStateQueue no more than once. 

Theorem 14 (Algorithm |2] terminates). Algorithm\^ terminates after the action of no more 
than Axn(n — 1)^/2 lowering operators 

Proof. In each of the n — 1 stages of Algorithm El the states that are added to currentState¬ 
Queue are LI of each other because of the conditions imposed in the algorithm. There are no 
more LI states in the given SU(n) irrep than the dimension of the irrep space. Thus, the 
total number of states that are added to currentQueue in each of the n — 1 stages is no more 
than Ax- No more that n{n — l)/2 lowering operators are applied on the states that enter 
currentQueue. Thus, each stage terminates after the application of Axn^n — l)/2 lowering 
operations. Furthermore, the algorithm terminates after n — 1 stages and the application of 
no more than AKn{n — 1)^/2 lowering operations. ■ 

Finally, we prove that the canonical-basis-states algorithm returns the correct output 
when it terminates. 


Theorem 15 (Algorithm El is correct). The SU{n) states 




yielded by AlgorithmlM are the canonical states of Definition^ 


(47) 


Proof. The theorem holds if the states yielded by Algorithm El have well dehned weights and 
have well dehned irrep labels. First we show that the weight of each state in the output of the 
algorithm is well dehned. Each state in the output is obtained either by enacting lowering 
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operators on the hws or by taking linear combinations of states that have the same weight. 
Linear combination of states with the same weights have well dehned weights themselves. 
Thus, all the output states have well dehned weights for SU(m) irreps for all 2 < m < n. 

We prove that the states have well dehned SU(m) irrep label separately for m = n and 
for 2 < m < n — 1. The correctness of the su(m) hws follows from Lemma [6l Every state 
in the output is a linear combination of states obtained by lowering from the constructed 
su(m) hws. Thus, every state is in the correct SU(n) irrep 

The algorithm (Line [8]) enacts raising operators on linear combinations of su(m +1) basis 
states at one weight until each of the raising operators annihilates the raised state. The 
su(m) state thus obtained are legitimate su(m) hws’s or possibly linear combinations of 
su(m) hws’s by construction. The uniqueness of the hws is guaranteed by the existence of 
the canonical basis^. Each of the canonical basis states is obtained by lowering from these 
su(m) hws’s using su(m) lowering operators and thus have well dehned irrep labels for all 
2 < m < n — 1. 

We have shown that the states yielded by the algorithm have well dehned values of irrep 
labels for su(m) algebras for all {m : 2 < m < n} and of su(m) weights for all 
{m : 2 < m < n}. Thus, these states are the canonical SU(n) basis states. This completes 
the proof of correctness of Algorithm [2l ■ 

We have proved that Algorithm [2] terminates and returns the canonical basis states on 
termination. The states constructed by the canonical-basis-states algorithm are employed to 
compute arbitrary SU(n) P-functions using an algorithm presented in the next subsection. 


C. P-function algorithm 

Our task is to construct the D-function 


^A(") ,...,A(3), A(2) ;A'(") ,...,A'(3), A'(2) 

for given labels of the SU(?7,) element fo(0) given by the 

parametrization O. The "D-function fl48|) is computed as the inner product between the state 


^A(") ,...,A(3),A(2) / 


(49) 
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Algorithm 3 "D-fnnction Algorithm 

Input: 

• n G Z"*" t> Algorithm constructs P-functions of SU(n) elements. 

• ff = {cji, W 2 , ■ ■ ■, w„ 2 _i} G > Parametrization of SU(n) transformation. 

• , ■ ■ ■, and ,..., i> Row Label. 

• ..., and A'^"'),..., A'(^) t> Column Label. 

Output: 


1 : 

2 : 

3: 

4: 

5: 

6 : 

7: 

8 : 

9: 

10 : 

11 : 

12 : 


procedure D(n, Sd, ,..., , A'(«),..., K '(^), A^”),..., A^^), A'(”),..., A'(2)) 

Construct V G GL{n, C) from parametrization 
if A(^) = A'(^) then 

'^a(")’ '’a(3)’aS') ^ using CANONICALBASISSTATES(n, A(n)). 


Ip 




A'(") ,...,A'(3), A'(2) 

Construct (’^( 3 )/a( 2 ) 


using CANONICALBASISSTATES(n, K'{n)). 
from 


^A(") ,...,a(3), a(2 ) / complex conjugation. 


Construct V{VL) 


^p 


Return P = ’ .^aO), a( 2 ) 

else 

Return T> = Q 

end if 

end procedure 


A'(") ,...,A'(3), A'(2) 

R(L!) 


A 


using ^ Y.j 

A'("b...,A'(3),A''(2)' 


Ip 


A'M ,...,A'(3)^A'(2) 


of Eq. flSTl) and the transformed state 


1 /( 0 ) 




A'(") ,...,A'(3), A'(2) /• 


(50) 


Algorithm [3] constrncts the fnndamental representation, i.e., the n x n matrix, Vij of the 
SU(n) element 1/(0)^. Then, the expressions for the basis states 


, A("),...,A(3),A'(2) 
XA(n) ,...,A(3), A(2) 


\ , A'("),...,A'(3),A''(2)\ 

'*/a'(") ,...,A'(3), A'(2) / 


(51) 


corresponding to the given labels are compnted using the canonical-basis-states algorithm. 
The basis states thus obtained are expressed as summations over products of creation and 
annihilation operators. V (O) acts on the boson realization by transforming each boson 
independently according to 

A ^ = E (32) 

k 
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where are the matrix elements of the nxn representation of The algorithm 

transforms the second basis state of Eq. flSTl) nnder the action of V (ff) by replacing each of 
the the creation and annihilation operators of the state according to Eq. fl5^ . 

The "D-fnnction is evalnated as the inner prodnct using the commutation relations fl22|) 
or equivalently by using the Wick’s theorem^. The correctness and termination of Algo¬ 
rithm [3] follows directly from Theorems [14] and [151 which completes our algorithms for the 
computation of boson realizations of SU(n) states and of "D-functions. 


V. CONCLUSION 

In summary, we have devised an algorithm to compute expressions for boson realizations 
of the canonical basis states of SU(n) irreps. Boson realizations are ideally suited for an¬ 
alyzing the physics of single photons, providing a tractable interpretation to basis states 
as multiphoton states, and to transformations on these states as optical transformations. 
Furthermore, we have devised an algorithm to compute expressions for SU(n) D-functions 
in terms of elements of the fundamental representation of the group. Our algorithm offers 
signihcant advantage over competing algorithms to construct "D-functions. Furthermore, our 
P-function algorithm lays the groundwork for generalizing the analysis of optical interfer¬ 
ometry beyond the three-photon leveU^^— . 

This work is the hrst known application of graph-theoretic algorithms to SU(n) represen¬ 
tation theory. We overcome the problem of SU(n) weight multiplicity greater than unity by 
modifying the breadth-hrst graph-search algorithm. Our procedure for generating a basis 
set can be extended to subgroups of SU(n). In particular, the boson realization of the hws 
of 0(2/c) and 0{2k + 1) irreps can be constructed along the lines of Lemma Graph- 

search algorithms can be employed to construct 0(n) basis states and "D-functions if the 
problem of labelling 0(n) basis states can be overcome. Our approach opens the possibil¬ 
ity of exploiting the diverse graph-algorithms toolkit for solving problems in representation 
theory of Lie groups. 
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Appendix A: Choice of subalgebra chain 

Our algorithms construct canonical basis states that reduce the subalgebra chain flT 6 |) . 
Other su(n) D su(n — 1) D • ■ ■ D su(2) subalgebra chains are possible and our algorithm can 
be generalized to construct canonical basis states that reduce other chains, as we discuss in 
this appendix. 

Each su(m) subalgebra of su(n), m < n is specihed by the sets of raising, lowering and 
Cartan operators that generate it. For a given sequence 

= (Ai) 

of m increasing integers, we can dehne the corresponding set of raising, lowering and Cartan 
operators 


*{ ,22 5 ^il ,23 5 * * * 5 5 ^ 22,23 ; • 

r^- • 


(Raising) 

(A 2 ) 

*{ 5 ^^3)^1 5 • • • 5 ^imAl ’ ^^3)^2 5 * 

C- ■ 

5***5 1 } 

(Lowering) 

(A3) 

{C^ 22,22 ^^3)^3 ^ 22 , 22 ; • • 

r ■ 

* ? I'^m 


(Cartan) 

(A4) 

that generate the algebra. Thus, each su(n) D 

n 

T- 1 

1 

• D su(2) subalgebra 

chain is 


uniquely specihed by the ordered sequences : m < n of integers, where 

/'”-!) = 4”-'’. ■ ■ ■. 4”F} c { 1 , 2 ,..., 

/("-2) = {4"-=>,4”-">,...,4”-4’} c 

(AS) 


/{ 2 ) 



4^^} C 
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Consider the example of su(2) snbalgebras of 5u(3). The three snbsets 


{^*1,2, 0*2,1, 0*1,1 — 0 * 2 , 2 } 

(A6) 

{0*1,3, 0*g,i, 0*1,1 — 0*g,g} 

(A7) 

{0*2,3, 0*g,2, 0*2,2 — 0*g,g} 

(AS) 


of the generators {Cjj ■ i,j E {1, 2, 3}} of su(3) generate three distinct su(2) algebras. Each 
of the three snbsets (1A6ll - (1A8|l can be labelled with a two-element snbset of the {1, 2, 3} and 
can be employed to dehne canonical basis states of SU(n). For instance, consider (A, k) = 
( 1 , 1 ) irrep of SU(3). The weight (A 2 , Ai) = (0, 0) is associated with a two-dimensional space. 
We can identify two basis states of this space by specifying the following: 

1. choice of su(2) algebra. For instance = {1,2}, which corresponds to the algebra 
generated by {Ci, 2 , <^ 2 , 1 , Ci,i - (^ 2 , 2 }, 

2 . sui, 2 , 3 ( 3 ) irreps label; = ( 1 , 1 ), su(i, 2 )( 2 ) irreps label: = ( 0 ) and ( 1 ) for the 

two basis states. 

3. sui, 2 , 3 (3) weights: ( 0 , 0 ), su(i, 2 )( 2 ) weights: ( 0 ). 

Another basis set of the A = (0, 0) space of su(3) irrep K = (1,1) is specihed by choosing a 
different su( 2 ) snbalgebra as follows. 

1. choice of su(2) algebra. For instance = {1,3}, which corresponds to the algebra 
generated by {Ci^g, Ci,i - Cg^g}, 

2 - sui^ 2 , 3 ( 3 ) irreps label; = ( 1 , 1 ), su(i^g)( 2 ) irreps label: = ( 0 ) and ( 1 ) for the 

two basis states. 

3- sui^ 2 , 3 (3) weights: ( 0 , 0 ), sU(i^g)( 2 ) weights: ( 0 ). 

Thus, different choices of subalgebra chain give us different basis states. 

In the main text, we have chosen the subalgebra chain (IT^ . Our algorithms can be 
modihed to account for other choices of subalgbra chain by choosing a different set of lowering 
operators in the basis-set subroutine. Thus our algorithms can be used to construct states 
and "D-functions in any of the bases that reduce su{m) subalgebra chains. 
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Appendix B: Connection to Gelfand-Tsetlin basis 


In this appendix, we detail the mapping between our SU(n) basis states and the canonical 
Gelfand-Tsetlin (GT) basis. The GT basis identihes each SU(n) irrep with a sequence of n 
numbers 


■ ■ ■ ) ^n,n) 

mk,n > mk+i,n VI < fc < n - 1, (B2) 

where the hrst label in the subscript is the sequence index and the second label identihes 
the algebra. The carrier space of every 5u(m) subalgebra is composed of disjoint su(m — 1) 
carrier spaces 

(B3) 

that obey the betweenness condition 


'^k,n ^ ^ 'l^k+l,n- 

Thus, each su(n) basis state \M) can be labelled by the GT pattern 


(B4) 


|A/) = 


/ mi^AT 1712,N ■ ■ ■ 

mi,2 1712,2 

V "^1-1 




(B5) 




where 


(B6) 


mk,e > mk,n-i > mk+i,e, 1 < k < i < n. 

The canonical basis states are eigenstates of the Gartan operators {Hi} flT^ as detailed 
in the following lemma. 

Lemma 16 (Gonnection to Gelfand-Tsetlin basis^). The canonical basis states are con¬ 
nected to the GT basis according to 


^AW,...,AW,A(2) / 


mi,jq 777-2,AT . . . 171n,N 

mi,is[-i ... 


\ 


(B7) 
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Every state \M) in the GT-laheling scheme is a simultaneous eigenstate of all 5\x{n) Cartan 
operators, 

He\M) = Xf\M), (B8) 


with eigenvalues 

e ^ /i+i £-1 \ 

^ mk,£ - 2 I '^k,£+i + ruk^e-i j , I < i < N -1. (B9) 

k=l \fc=l k=l J 

Thus, the canonical basis states of Def. |4] is uniquely mapped to the GT basis. 

Furthermore, the weights A^ are also mapped via the boson realizations to differences in 
number of bosons at sites ^ and £ + 1. Hence, the difference 

z/£+i = ^ mk/ - 2 I X] mk,i+i + ^ mk,£-i \ (BIO) 

k=l Vfc=l k=l / 

in the number of bosons at sites ^ + 1 and ^ of the boson realization of a basis state is also 
connected to its GT pattern. Once we recall the total number of bosons in the system is 
+ 1^2 + ■ ■ ■ + = Nk, one can then invert the differences and recover Vp in term of the 

mk^i-i- Thus the canonical GT basis states are connected to our SU(n) basis states. 
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